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We investigate both numerically and analytically the behaviour of a spin-1 antiferromagnetic (AFM) 
isotropic Heisenberg chain in an external magnetic field. Extensive DMRG studies of chains up to 
N = 80 sites extend previous analyses and exhibit the well known phenomenon of the closure of 
the Haldane gap at a lower critical field H C1 . We obtain an estimate of the gap below H C1 . Above 
the lower critical field, when the correlation functions exhibit algebraic decay, we obtain the critical 
exponent as a function of the net magnetization as well as the magnetization curve up to the 
saturation (upper critical) field H C2 . We argue that, despite the fact that the SO(3) symmetry of 
the model is explicitly broken by the field, the Haldane phase of the model is still well described by 
an SO(3) nonlinear cr-model. A mean-field theory is developed for the latter and its predictions are 
compared with those of the numerical analysis and with the existing literature. 



I. INTRODUCTION. 

Spin chains, ladders and, more generally, low- 
dimensional quantum spin systems have become a sub- 
ject of considerable theoretical and (more recently) ex- 
perimental interest Q. This has been due in the recent 
past to the connection that 2D AFM Heisenberg models 
have with the Cu—O planes in the new ceramic supercon- 
ducting materials [^), but spin chains and ladders have 
become rapidly a subject of autonomous interest in their 
own. 

Back in the early 80's Haldane || put forward his by 
now famous "conjecture" according to which half-odd- 
integer AFM Heisenberg spin chains should have a gap- 
less spectrum, and hence be critical at T — with al- 
gebraic decay of correlations while integer spin chains 
should be gapped with exponential decay of correlations. 
The conjecture is based on the presence, in the effective 
action that, in the continuum limit, maps the chain onto 
a (1 + 1) nonlinear Sigma model [NLaM] H, of a topo- 
logical term (basically a Pontrjagin index Q-0|), that is 
effective only when the spin is half-odd-integer (for de- 
tails see, e.g. pi). It was later proved rigorously || that 
a (1 + 1) NLaM without topological terms has a unique 
ground state and is gapped with an excitation spectrum 
consisting of a triplet of massive bosons, as well as that 
the topological term makes the model gapless Q. This 
implies of course that all half-odd-integer spin chains fall 



into the same universality class of the spin-1/2 chain, 
whose exact solution is well established by the Bethe- 
Ansatz |lC| , |ll| ], and that is known to be gapless. Were it 
not for the presence of the topological term, the descrip- 
tion of half-odd-integer spin chains in terms of effective 
nonlinear a model [NLctM] actions would lead then to 
predictions in blatant disagreement with the results of 
the Lieb-Schultz-Mattis [LSM] theorem (l2), which pre- 
dicts a (non Goldstone) gapless sector in the excitation 
spectrum in the presence of a unique ground state. All 
this applies to systems with translationally- (and parity)- 
invariant couplings. Gapfull half-odd-integer spin chains 
can occur when these invariances are explicitly broken 

0. 

The conclusions of LSM do not apply to integer- 
spin chains, though, and there a unique ground state 
and a gap (the "Haldane gap") may well coexist. In- 
deed, neutron-scattering measurements on quasi one- 
dimensional compounds |l3f| have established the pres- 
ence of the Haldane gap in a class of integer-spin sys- 
tems. Also, a rigorous example of an S = 1 model with a 
gapped unique ground state has been exhibited (among 
other rigorous results) in a seminal paper by Affleck et 
al. @. 

In more recent years there has arisen a growing interest 
in the behaviour of spin chains and/or ladders in an exter- 
nal (uniform and/or staggered jl5|]) static magnetic field. 
Actually, the magnetization curve of the spin one-half 
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AFM Heisenberg chain had been studied already by Grif- 
fiths back in the early Sixties 16 1 . The main result of Ref. 
Jill was that the magnetization m per site, as a function 
of the external field H, grows steadily from H — (with 
a finite slope corresponding to a finite initial susceptibil- 
ity, as expected) to the saturation value of to = 5 = 1/2. 
The latter is reached at an upper critical field H C2 given 
by (absorbing the Bohr magneton and the spin g-factor in 
the definition of the field): H C2 = AS J = 2 J for S=l/2 
(J being the AFM exchange constant) with infinite slope: 
5 — to oc (H C2 — H) 1 / 2 . A continuous (in the thermody- 
namic limit) change of the magnetization as a function 
of the applied field implies of course that there are no 
gaps in the spectrum towards excitations from a value 
of the spin to a neighboring one. The 5 — 1/2 model 
remains therefore gapless from H = up to H = H C2 , 
where the ground state becomes fully polarized. In this 
connection, it should be mentioned that the opposite phe- 
nomenon, namely a field-induced gap has been observed 
Jl7f in neutron-scattering experiments on what can be 
considered as a (quasi) one-dimensional S = 1/2 spin 
system, namely Cu benzoate. It has however been ar- 
gued convincingly jl8j that this is possible only in the 
presence of various types of anisotropies, leading to an 
effective internal staggered field (see 1513] for more de- 
tails). 

As already stated, integer-spin (isotropic) AFM chains 
are in a gapped "Haldane" phase at zero field. This ap- 
pears to be true Jll| also for XXZ models, i.e. for more 
general model Hamiltonians of the form: 



N r i i 

H = J E { A ^ 5 m + 2 (^ +5 r+i + SiS+ +1 ) | (1 



for a chain of N sites (we don't specify here boundary 
conditions (see jl9| for details)), where A > is an 
anisotropy parameter in a range that includes the (ex- 
actly soluble, A = 0) XY model, the (isotropic, A = 1) 
Haldane phase as well as, for A > 1, uniaxial Ising-likc 
phases. 

Sticking for the time being to the isotropic, A = 1 case, 
and according to the analysis of ||, the lowest excited 
state is a triplet of massive, 5 = 1 bosons. In the absence 
of an external magnetic field the latter are degenerate and 
separated from the ground state by a finite energy gap 
Aq. The addition to the model of an interaction with an 
external field leads to new and interesting phenomena. 

As observed by Affleck [M , the Zeeman splitting due 
to the field will cause one of the branches of the spectrum 
to lower its energy. The gap will remain "robust" (and 
the magnetization will remain zero) up to a lower critical 
field: H Cl = Ao, when the lower branch of the spectrum 
will cross the ground state. At this point, a "Bose con- 
densation of magnons" [ElJ should take place, and the 
spectrum should become gapless. This implies of course 
that the magnetization should remain zero for fields up 



to the lower critical field, while it should rise to the satu- 
ration value at some upper critical field H C2 > H Cl . The 
value H = H C1 of the field is then the one that leads 
to the closure of the Haldane gap in the excitation spec- 
trum. Parenthetically, it should be remembered that the 
existence, for integer spins, of a lower critical field lead- 
ing to the closure of the Haldane gap had already been 
predicted previously by Schulz p2] j. 

The Lieb-Schultz-Mattis Theorem JE|, originally for- 
mulated in zero external field, has been generalized to 
both spin chains [^3|,^l| and ladders (24j| in the presence 
of a field. The main results of (2l]] for spin chains (but 
see also p3] for a generalization) are that, in the presence 
of a field: 

i) If translational symmetry is not broken, then the 
ground state for a spin-5 chain is gapless unless the mag- 
netization per site to obeys the "quantization condition" : 



5 — m G Z, an integer 



(2) 



This implies that gaps in the spectrum, and hence 
plateaus in the magnetization curve can (but in principle 
need not, of course) be present when the above "integer 
quantization" condition is met. Moreover: 

ii) "Fractional quantization" (completing an overall 
scenario that resembles closely |2l]] that of the (integer 
and fractional) Quantum Hall Effect), i.e.: 5 — to = half- 
odd-integer, may lead to additional "fractional" plateaus 
accompanied, however, by spontaneous breaking of the 
translational symmetry of the ground state. 

Although plateaus (both integer and fractional) are 
not forbidden at quantized values of the magnetization, 
whether or not they do actually occur depends very much 
on the details of the model. For 5=1 the only admissible 
gaps are at to = (with a width: H Cl w A , the Haldane 
gap) and (trivially) at m = S = 1 for H > H C2 , corre- 
sponding to the fully polarized state. In order to have 
gaps at intermediate values one needs therefore 5 > 3/2. 
The existence of a gap at to = 1/2 has been proved 
by Tasaki (2^] and numerically found [§1 26 only in the 



presence of a strong single-ion anisotropy. Some numeri- 
cal |24[ and experimental |2?j evidence of magnetization 
plateaus has been reported for 5 = 1/2 (even- legged) 
spin ladders. 

We will restrict ourselves from now on to the case of 
5=1 AFM isotropic Heisenberg chains. For these sys- 
tems the magnetization is supposed to rise continuosly for 



< H < H C2 from zero to saturation. 



The model be- 



ing gapless in this range of fields, the transverse spin-spin 
correlation function (e.g.: (5q5^) if the external field H 
is in the z-direction) is expected to decay algebraically: 
(5q5^) oc (— l)'Y~ r ' with some critical exponent to Af- 
fleck Q has suggested that the ground state above H Cl 
should be regarded as a Bose condensate of low-energy 
bosons and has proposed an effective Ginzburg-Landau- 
type field theory to describe the system near the transi- 
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tion. His predictions are that r\ should exhibit the uni- 
versal behaviour: r\ ~ 1/2 near H Cl and that the induced 
magnetization should vanish in a nonanalytic way at H Cl , 
namely: m oc tJH — H C1 for H > ff ci . This prediction 
seems to be reasonably well supported by numerical find- 
ings |2£| . In the same reference it was also inferred from 
the numerical data that r\ w 1/2 as well when ff — > i/ C2 . 
It is unclear whether such a prediction can be actually 
extracted firmly from the available numerical data. 

The motivations of the present paper are twofold: we 
wanted to extend and to put on firmer bases the existing 
numerical analyses of the model, which has proven to 
be possible by a systematic use of the DMRG that has 
allowed us to analyze samples with an higher number of 
sites (up to N = 80) than down hitherto, and to try and 
understand the Haldane phase in terms of an appropriate 
effective field theory. 

The plan of the paper is as follows: In Sect. II we give 
the details of the numerical procedure followed in the 
study (at T = 0) of finite samples and report the results 
for the Haldane gap, the critical fields, the magnetiza- 
tion curve and the critical exponent for the (transverse) 
spin-spin correlation function. In Sect. Ill we discuss, at 
the mean field level, the predictions of an effective field 
theory based on the 0(3) nonlinear a model (NLcrM) for 
the Haldane phase of the chain, comparing the predic- 
tions of the theory with the numerical results. Sect. IV 
is devoted to a brief discussion and to the conclusions. 

II. NUMERICAL RESULTS 

The Hamiltonian we consider here is: 

n = j^s.-s^-tf^s? (3) 

(ij) ■>■ 

where: S? = S(S + 1), S = 1 (we set h = 1 from now 
on), the sums run over a chain of N = L/a spins, L 
being the length of the chain and a the lattice spacing. J 
is the AFM coupling between nearest-neighboring spins 
and H (which includes factors such as Bohr magneton, 
gyromagnetic factor and the like) is a static and uniform 
external magnetic field pointing in the z-direction. Here 
and in the rest of the paper we shall assume periodic 
boundary conditions (PBC's). 

We have employed the standard DMRG algorithm 
p9| [j2| to calculate the magnetization curve and the 
static transverse spin-spin correlation function: 

G xx {r) = (S X S%) (4) 

(r = j, in units of a with j an integer, < j < N). From 
the latter one can extract both the value of the correla- 
tion length (and hence of the spin gap) in the Haldane 
phase and, as a function of, e.g., the magnetization, the 



critical exponent governing the algebraic decay of corre- 
lations in the gapless phase. We have adopted a finite- 
chain "zip" algorithm which increases considerably the 
precision of the resulting data, though of course at the 
price of an increase in computer time. In the whole of 
this Section we measure energies in units of J, i.e. we set 
J = 1. 
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FIG. 1. a) The hystogram of the magnetization curve for 
N = 32 spins, b) Continuous interpolation of the data in Fig. 
1-a. See text for details. 

The hystogram of the magnetization for N = 32 spins 
is reported in Fig. 1-a, while Fig. 1-b represents a con- 
tinuous interpolation of the data. The curve in Fig. 1- 
b passes through the midpoints of the plateaus of the 
hystogram and is an indication of what the magnetiza- 
tion curve should look like in the thermodynamic limit 
N -> oo. 

The magnetization remains zero up to a lower criti- 
cal field of H Cl = 0.412232 (in our units). Identifying 
H Cl with Ao, Ao being the spin gap, this is a close up- 
per limit to the N — > oo extrapolations [^9| that give: 
Ao/ J—0. 40450(2). Similar results have been reported 
also by Golinelli et al. |53|. Above the lower critical field 
the magnetization increases steadily until it saturates at 
the upper critical field H C2 /J — 4 (equal to 45 for S = 1). 
From the data of Fig.l, it emerges that the numerical 
analysis reproduces in a completely nonambiguous way 
the singular square-root behaviour of the magnetization 
near the upper critical field. The data are also consistent 
with a similar square-root vanishing of the magne- 
tization when the lower critical field is approached from 
above, but seem to indicate that the low-field singular be- 
haviour sets in in a much narrower range of fields than the 
one near H C2 . The values of (— l) r G xx (r) are plotted in 
Fig. 2 for a chain of N = 64 spins and H < H C1 . PBC's 
make of course the curve symmetric around N/2 = 32. 
The data are well fitted, taking PBC's into account, by 
a curve of the form: 

,- irC - W ^ {_ + __} (5) 

with: 

A = 1.07508, £ = 6.08304 (6) 

both independent of the field up to H = H Cl . The 
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form of G xx (r) reported here corresponds to the (sym- 
metrized) asymptotic behaviour of the modified Bessel 
function JjiJ Kq(t/£), which is the form of the spin-spin 
static correlation function predicted by the continuum 
theory to be discussed in the next Section. Also, defin- 
ing the spin- wave velocity c as: c — Ao£ (Ao£/H if we 
reinstate Planck's constant into its proper place) we find 
c ~ 2.4 in units of J a (or Ja/K). The continuum the- 
ory (see below) predicts instead c — 2 (for 5 = 1 and 
in the same units). The discrepancy may be attributed 
both to the approximations involved in performing the 
continuum limit and to finite-size effects that may affect 
the numerical estimates. 
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FIG. 2. A plot of (-l) r G xx {r) (r = j, in units of a, where 
j is an integer) for a chain of N — 64 spins and H < H C1 . 
The points are the numerical data and the continuous curve 
is the interpolation with the analytic form of Eq.(5). 

For H Cl < H < H C2 the correlation function is ex- 
pected to exhibit an asymptotic algebraic decay, i.e. to 
behave asymptotically as: 



(7) 



Here too the second term accounts for the PBC's and 
i] = r)(m). According to Moreo |35) and to Sakai and 
Takahashi |^8| the most efficient way to calculate r\ goes 
through the evaluation of the (staggered) structure fac- 
tor: 



Sjy — 



N-l 

E 

r=l 



exp{ikr}G xx (r)\ k=7T 



(8) 



A power-law behaviour as in Eq.(Q) implies: Sn ~ 
N 1 ^ v for large N, and one obtains the critical exponent 
as: 



7y(m) = 1 — 



ln(S N /S M ) 
In (N/M) 



with the most reliable values being obtained for the high- 
est M, N available. 

In Fig. 3 we report our numerical results for ry together 
with some interpolations. The data have been obtained 
by employing about eight runs for each value always using 
the "zip " algorithm twice and keeping between 190 and 
234 states in the reduced density matrix. The longest 
chains studied here have up to N — 80 spins. 
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FIG. 3. The crtitical exponent r\ (see Eq.(7)) as a function 
of the magnetization for chains up to TV = 80 spins. The 
interpolating curves should be taken mainly as guides for the 
eye. 

For unOwe obtain extrapolated values of r\ slightly 
above the theoretical prediction [^0| of r\ = 1/2. This 
may be due to uncertainties in the interpolation scheme 
and again to finite-size effects. For m — ► 1, i.e. near 
saturation, our data seem to be consistent with a finite 
limiting value of r\ at saturation. 

Our numerical results indicate that the system remains 
in the Haldane phase for fields up to H Cl , while it be- 
comes (quantum) critical for higher fields when a nonzero 
magnetization develops. In the next Section we will try 
to set up a consistent field-theoretical description of the 
Haldane phase. 



III. EFFECTIVE FIELD THEORY FOR THE 
HALDANE PHASE 

That AFM Heisenberg spin chains can be mapped, in 
the continuum limit, into 0(3) NLu models was shown 
for the first time by Haldane || and such "Haldane maps" 
have been widely used since. As the procedure is by now 
well known, we summarize here only its main points. 

In order to use the partition function as a generating 
functional for (Euclidean) Greeen functions, we will re- 
place for the time being the Hamiltonian of Eq. (|3|) with: 



4 



H = J^Si-Sj-^Bi-Si, J>0 (10) 

(y> « 

where Bi is an external site- and possibly time-dependent 
magnetic field. Here too PBC's will be assumed through- 
out. 

The (canonical) partition function can be written in a 
standard manner as a path integral over spin coherent 
states as: 



Z = J [Vfl\5(fl 2 - 1) exp-t-Sfi} 

where the (Euclidean) action is given by: 

P 

S E = S wz + J H[Sl(T)]dT 



(11) 



(12) 



f2 is a classical vector field of unit magnitude, 
W[0(t)] = JS 2 toiir) ■ n j( T ) -^B s (t). ^(r) 



(13) 



and Swz is the Wess-Zumino part of the action given by: 



Swz = is 



rP 

/ dr[l — cos 8, 
„ Jo 



ir)]Mr) (14) 



where 0i and <pi are the polar coordinates of the 
parametrization of tti. Following Haldane, we decom- 
pose r2j into a staggered Neel field and a fluctuation 
field L as: 



a = (-i)Wi-« J i? + <«i< 



(15) 



and the constraint fif = 1 is enforced by: nf = 1 and 
rij • I, = 0. 

Performing then the continuum limit followed by a gra- 
dient expansion and integrating out the fluctuation field, 
one ends up with the following expression for the parti- 
tion function expressed as a path integral over the Neel 
field alone: 



Z = J [£>n]<5[n 2 - 1] cxp[-S £ 



with the Euclidean action: 

l p 



Se = j ' dx J dr£E(x,T) 



(16) 



(17) 



o o 



and 



Ce(x,t) = 1 [\\drn\ 2 + C |d,n| 2 } — (B 2 - (B • n) 2 



where g = 2/S and c = 2JSa are the bare coupling 
and spin-wave velocity, respectively. This form of the 
Euclidean Lagrangian has already been derived by other 
authors p6| -f38|. The last term in £ B (x, r) is basically 
geometric in nature, being linear in time derivatives, and 
comes from the variation of the Wess-Zumino part of the 
action. The Berry phase instead, that comes also from 
the Wess-Zumino action, can be safely neglected as we 
are dealing here with integer spins. 

An additional term that is present, e.g., in the deriva- 
tion of the effective action in Ref. |Q gives rise, in the 
continuum limit, to an integrated total (space) deriva- 
tive that, as a consequence of the assumed PBC's, does 
not contribute to the overall Euclidean action. However, 
it can affect in a crucial way the evaluation of the spin 
Green functions. Indeed, using Eq.(|l5|), the spin Green 
functions are given in the exact (lattice) version of the 
theory by M: 



— B • nx<9 r n 



G 



ap 



(r) = (t t {S?(t)SP (0)}) 



5 2 (-l)^(T T {<(r)nf(0)} 
S 2 a 2 (T T {Jf(r)Zf(0)} 



(19) 



to leading order in 1 and with the expectation values of 
the cross terms vanishing by symmetry. According to 
the analysis of the first of Refs. ^] the second term 
in Eq.(|lS|) (a two-magnon contribution) costitutes a mi- 
nor correction to the dominant, staggered part and, in 
the static (i.e. equal-time) limit it decays asymptotically 
faster than the staggered part. So, we will approximate 
the spin Green functions as: 



Gtf(r) ~ S 2 (-l)^ (T r {n?(r)nj(0)}) (20) 



having in mind that, due to the staggering prefactor, the 
small-q behaviour in Fourier space of the Green functions 
for the n-field will correspond to that of the spin Green 
functions near q = n. 

One may wish also, for reasons of internal consistency, 
to derive the Green functions (actually the truncated 
ones) directly from the partition function via functional 
differentiation, i.e. as: 



S 2 Z 



8B?(t)8BUQ) 



(21) 



(18) 



This can be done safely, however, only before the contin- 
uum limit is taken, as can be seen easily from a careful 
analysis of the Zeeman term in the Hamiltonian (p^). 
With the substitution ( fl5| ) the latter becomes: 

S^B ; (t) • n,( T ) = ^{(-IWt) • n,(r) 

i i 

+ aBi(r) • li(r)} (22) 
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and it is pretty obvious that the double functional differ- 



entiation of Eg. (J2l|) will reproduce the r.h.s. of Eq.(19). 
The first term in Eq. (|2^) is however precisely the one that 
maps, in the continuum limit, into the already-mentioned 
irrelevant (in that limit) total derivative term. So, if func- 
tional differentiations were performed (erroneously) after 
taking the continuum limit they would only reproduce 
the nonstaggered part (with the site indices i,j replaced 
by continuos variables, of course) of the Green functions 
(|l9|). This we have checked also by direct calculation, 
performing the functional differentiation both before and 
after having done the Gaussian integration over the fluc- 
tuation field. 

Returning now to the evaluation (in the continuum 
limit now) of the partition function, the NLerM constraint 
can be implemented by re-expressing the partition func- 
tion as: 



Z = J [Vn] 



2^ 



CXpI-.Se//] 



where now: 



L 



S eff = J ' dx J dT £eff{x,T) 



(23) 



(24) 



o o 



and 



B a (x,r) = -^^exp[i( qa ;-n n T)]B Q (q,n) (28) 

L 

B a {q,n) = J dx J dTexp[-i(qx-n n T)}B a (x,T) (29) 



and similarly for the components of n and for the field 
A. The frequencies here are Matsubara-Bose frequencies: 
fl n — 2nn/(3. Note that reality of B, n and A implies: 



as well as: 



and 



B a (-q,-n) = {B a (q,n)r 



l a (-q,-n) = {n a (q,n)y 



X(-q,-n) = {X(q,n)Y 
Some long algebra leads then to: 



(30) 



(31) 



(32) 



S ?ff = E E ""fa- ny(K- l T~<{q, n; q , n , )rP(q', n') 

q,q' n,n' 

,2 



-^EI^)l 2 +*°) 

q , n , a 



(33) 



where 



C eff (x,T) = C e (x,t) - i\{x, T ){n 2 {x,T) - I) (25) 

with A(x, t) an auxiliary (real) Lagrange multiplier field. 
The total effective action is now quadratic in the field n 
and is given by: 

L L 

S e ff — J dx J dr J dx' J dr' 

oooo 

{n a (x, T)(K- 1 ) a '> \xt- a/rOnV, r')} 

L 

+ j dx J dr ^~^-B 2 {x,t) +a(x,r)| (26) 
o o 

where: 

{K-^txT-x'T 1 ) = 

{-hU* + cdl]5^ + -L B"(x,t)B^(x't') 
I 2g c 2gc 

+— e a ^B s (x,T)d T - i\(x,T)5 a A 6(x - x')6(t - t') (27) 

gc J 

and we can proceed to integrating out n, leaving a result- 
ing effective action depending only on B and the auxil- 
iary field A. The integration involves the evaluation of 
the functional determinant of K 7 and this is most easily 
done in Fourier space. 

Fourier transforms will be defined according to: 



{K- X ) a ~i(q, n; q', n') = <W S qq 'S^^[n 2 n + c 2 q 2 ' 



2gc L 



1 1 
1 1 



V B a (k + q, m + n)B~<{k + q' , m + n') 



gc (/3L)2 
1 



-i5 an >- 



n n> e a ^B\q-q',n-n') 
7 X(q- q',n- n') 



(34) 



It is clear that the Fourier transforms of the (nonlocal 
in the presence of a general external field and/or La- 
grange multiplier A) Green functions of the n field are 
given by: 



(n a (q, n)n!{q', n')) = K ai {q, n; q', n') 



(35) 



We can perform now the Gaussian integration over the 
field n with the result: 



Z = J [£>A]exp{-S[A;B]} 



(36) 



where 



5[A;B] = Trln^" 1 ) 



1 1 

2g~cJZ 



]T \B Q (q,n)\ 2 +*A(0,0); 

q,n,a 

(37) 
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Trln^- 1 )^ Y.iHK-^iq^n-q.n) (38) 

From now on we will limit ourselves to the case of a 
static and uniform magnetic field, i.e. we will set: 



B a (q,n) = pL5 qfi 5 nfi H° 



(39) 



In that case: 



{K-^^n-q^n') = (K^(q,n;q' ,n') 

f 

0L) 



- iS ar f 7 ^ F ^\(q-^,n-n') (40) 



where K& is the diagonal (in Fourier space) kernel given 
by: 

Kd 1 (.Q' n ''<l'' n ') = J^Snn'Sqq'G^iq^n) (41) 

with: 

+ —H a H"< + —Q, n e al5 H s (42) 
2gc gc 

Looking for the saddle point (s) of the effective action, 
we find that there exist static and uniform ones of the 
form: 



X(q,n) = X/3LS q fiS n fi 



(43) 



implying: 



(K-^fen^V) = ^<WV0 _1 (9>») ( 44 ) 



and: 



(g-'THq, n) = {Go^iq, n) - iX5 a ~> (45) 



At the saddle point the Green functions are diagonal (in 
Fourier space) and given by: 

(n a (q,n)n^(-q,-n))=g^(q,n) (46) 

The saddle-point equation reads then: 

where "Sp" stands for a trace over the vector indices only. 
Frequency summations of the form: 



(48) 



are performed according to a standard recipe as: 

^F(Sl n ) = ^F{z) P {z)^ (49) 



where the contour T goes from — oo — i0 + to +oo — i0 + 
and backwards from +oo + i0 + to — oo + i0 + and: 



p(z) = {cxp[if3z} - 1} 1 



(50) 



is a meromorphic function having simple poles at: 
z = 0„ = 27m//3 with residue Closing the first 

integral in the lower half-plane and the second one in the 
upper half-plane the resulting integrals can be evaluated 
with the aid of the residue theorem. 

Let's begin by examining the zero-field case: H a = 0. 
Then: 

{Q- l r{q, n) = { ^{fin + C V} iX ] ( 51 ) 



and the saddle-point equation becomes: 



1 x - 1 



65c 



r 



-2iXg 



(52) 



Under the assumption that £ -2 be positive (i.e. that 
the saddle point be actually purely imaginary and in the 
upper half-plane), the integrand has simple poles at z = 
±icy* 'q 2 + £~ 2 . Hence: 



-r 



6gc 



3.9 



cothj^cv^ + F 5 } (53) 



In the (thermodynamic) limit L — > 00 the momentum 
integral requires the introduction of an ultraviolet cutoff 
in momentum space and we obtain: 







and, in the zero-temperature limit (3 — > 00: 

A 

2tt- 





(3g/2n — 3/ttS in our case) leading to: 



£0 = £(ff = 0) = - sinh(2^/3. 9 ) 



and to the diagonal propagator: 
ai (q,n) = 
G(q,n) = 



g^(q,n)=5^g(q,n) ] 
2gc 



(55) 
(56) 

(57) 
(58) 



The ultraviolet cutoff will be taken here as a free parame- 
ter to be determined by fitting the results to the existing 
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data for the spin gap. At the saddle point the poles of 
the Green functions (that are all equal and diagonal) are 
a degenerate triplet of massive excitations with energy: 



and a zero- field gap: 

A = c^ 1 w A exp{ -ttS/3} 



(59) 



(60) 



in the large S limit. This is in agreement with the known 
exact results on the model [|| and with established esti- 
mates of the behaviour of the gap for large S Ex- 
plicitly, performing the frequency summation, the static 
propagator is given by: 



£(9) = ^£S(<7,«) = 



(61) 



and, in real space: 

0(x) = J ^ exp{iqx}G(q) = 



K (\x\/£ ) « exp{-|z|/£o}/VM (62) 



for \x\ ^> in agreement with the results of Sect. II 
and with previous predictions p0| . Therefore £o can be 
identified with the zero-field correlation length. 

In the general case, fixing the direction of the field 
along the z-axis: H a — H5 a3 , the matrix of the inverse 
propagators has the form: 



(63) 




where: 



A = A(q,Q n ) 
B = B(q,fl n ) 
C = C(q,n n ) 



2gc 

* H ' 2 
A H 

2gc 
— Hfl n 



2 2l 

C q } 



(64) 
(65) 
(66) 



Hence: 



Q = 



A 


c 


\ 


A 2 +C 2 


A 2 + C 2 


C 


A 





A 2 +C 2 


A 2 +C 2 








* I 



(67) 



and: 



Sp(G) 



2A 



A 2 + C 2 



1 

B 



(68) 



Proceeding as before in the evaluation of the frequency 
sums, we find: 



y/q 2 + S,- 2 + (H/c) 2 I 2 



2A 



and: 

-y 

3^ A 2 + C 2 

n 
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y/q 2 + r 2 + (H/c) 2 



coth 



^q 2 + ^ 2 + (H/c) 2 + H) 



coth 



^q 2 +i- 2 + {H/c) 2 -H) 



(70) 



For H — > we recover the previous case. In the zero- 
temperature limit the self-consistency equation reduces 
to: 



1 



(71) 



Letting L — > oo and introducing the same ultraviolet 
cutoff A as before wc find: 



1=^ 
2tt 



dry 



(72) 



Comparison with Eq.(|55|) tells us immediately that the 
value of £ -2 for H ^ is connected with that at H = 
by: 

r 2 (#) + (ff/c) 2 =£o- 2 (73) 

which means that £ _1 will vanish at a lower critical field 
H Cl given by: 



H C1 = c& 1 = A 



(74) 



Therefore i? Cl will be equal to the zero-field gap, and the 
divergence of £ for H — > iJ cl will be the signal of a phase 
transition at H = H Cl .Explicitly: 



(75) 



Notice, however, that at T > £ can never diverge, 
since for £ -2 = one of the hyperbolic cotangents would 
develop an infrared singularity which is not present at 
T = (strictly). 

We turn now to the Green functions. When evalu- 
ated at the saddle point, besides the diagonal ones, there 
appear also Green functions correlating the x and y com- 
ponents of the NLaM field n. We have: 
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(n x (q,n)n x (-q,-n)) = (n y {q,n)n v (-q, -n)) = 
(n x (q,n)n v (—q, —n)) = —(n v (q,n)n x (—q,—n)) ■ 
while: 

(n z (q,n)n z (-q,-n)) 



A 



A 2 



C 2 
C 



S = TrlniR- 1 } - — /3LA 2 
25c 



A 2 + C 2 



-and: 



1 

B 



(76) 



Interchange of the x and y directions in the internal 
space can be achieved by a rotation of 7r around, say, the 
x = y axis. This will result in a symmetry only if one 
sends simultaneously h — ► —h, and this explains the way 
the off-diagonal Green functions coupling n x and n v are 
related and can be deduced from one another. 

Notice that, being odd in frequency, the static 
(i.e., equal-time) off-diagonal correlation functions: 
(n x (q,T) n y(-q,T)) = 1//? n)n»(-g, -n)) = 
— (n v (q,T)n x (—q,T)) will vanish identically, as they 
should. There are therefore dynamical but not static 
off-diagonal correlations. 

Turning now to the elementary excitation spectrum, 
we observe that the degenerate triplet of excitations at 
H = gets splitted into a longitudinal branch with en- 
ergy e = eo(q) and two transverse branches with energies: 



4^ 



+ c 2 q 2 ± H 



(77) 



and gaps A = c^ 1 and A± = c^ 1 ± H, which seems 
to be appropriate ptj| for a triplet of S — 1 magnons in 
an external field. The gap of the lower branch will van- 
ish linearly as H approaches H Cl . Thus one mode will 
go "soft" at H Cl , and this is again a signal of instability. 
This result seems to be in complete agreement with the 
picture put forward some time ago by Affleck |2Q] accord- 
ing to which the transition at H Cl is due to a ID Bosc 
condensation of (soft) magnons. As it stands, however, 
the present mean-field theory cannot be extended in a 
straightforward way beyond the lower critical field H C1 . 
The static correlation functions (to be compared with 
those calculated numerically in Sect. II) can be derived 
immediately from the Green functions. In particular, for 
the n x — n x correlation function we obtain: 

(n x (q,T)n x {-q,T)) = 



-to 2 
(78) 



i.e. the static correlation function will be the same as 
in the zero-field limit, with the same correlation length. 
Therefore, at H ^ the physical gaps and the "gaps" in 
the static correlation function will have to be treated as 
related but distinct objects, coinciding only in the zero- 
field limit. 

We turn now to the calculation of the magnetization. 
The Gibbs free energy G is given by: G — (3~ 1 S where, 
at the saddle point: 



A^ = 



2,-2 ,_ H 2 



(79) 



(80) 



Notice that (cfr.Eq.©) at T = 0: A 2 = A 2 and is field- 
independent, but it need not be so at finite temperatures. 
Explicitly, and again at the saddle point: 

{K-^lXqn^n') = ^-W-W^ 1 )^ n) (81) 



Or 1 )^ (<?, n) = am— {a 2 + cV + r 2 )} 

— 6 a ' f S a ' 3 H 2 + — 0,e Q73 
25c 



1 



gc 

„2 „2 1 A 2 



5«-t-L- }tf + c 2 q 2 + A 2 - H 2 (l - 6 a3 )\ 
2gc 



H , 
— fine" 73 



(82) 



M 



dH 
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It follows then that the total magnetization is given by: 

(83) 

$°n QA 2 

(84) 



Now: 

d(Q- 1 )°n 
dH 



A a7 Q 



gc 



gc 



It turns out therefore that the overall coefficient multi- 
plying Qjjj- vanishes identically at the saddle point. The 
average magnetization per site: m = M/L is then given 
by: 

2H ^n 2 Jgc~A 



(3Lgc^ A 2 +C 2 

1 H Q2_ c 2 (?2 + r 2 ) 



PL (gc) 2 



q , n 



AH ^ 1 ^ 



A 2 + C 2 



Z+Z- 



(85) 



where z± = ie±(q) and we have used c 2 (q 2 + £ 2 ) 
— z+z-. The frequency sum is easily evaluated as: 



R ^ 







2H r th 2^ 



y/c*q 2 + A 2 -H 



- coth \ ~/3 \yjc 2 q 2 + A 2 + H 



(86) 
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Hence we obtain the cutoff-independent result: 



= {coth[i/3 \Vc 2 q 2 + A 2 - H 



■ coth <j -13 



y/c 2 q 2 + A 2 + H 



(87) 



which has to be considered together with Eq.([47|) and 
Eqs.(|68|)-([70|) that determine £ and hence A as a function 
of field and temperature. As a function of H, m is odd, 
nondecreasing and vanishing linearly for H — > as long 
as T > 0. However, the two terms in curly brackets cancel 
exponentially against each other when (3 — > oo, and we 
obtain the result that, as long as: H < H ei holds, i.e. 
in the whole Haldane phase < H < H C1 and at T = 
the magnetization vanishes identically, as it should. This 
completes the discussion of the mean-field theory for the 
0(3) NLtrM in the Haldane phase. 



IV. DISCUSSION AND CONCLUSIONS 

In Sect. II we have presented numerical results for 
finite chains. Our results show in a nonambiguous way 
the persistence of a Haldane (gapped) phase up to a lower 
critical field that can be identified with the spin gap as 
well as critical behaviour, with algebraic decay of correla- 
tions, for higher fields up to a saturation field. The shape 
of the magnetization curve and the values obtained there 
for the spin gap and the critical exponent above the lower 
critical field are in good agreement with previous work, 
both numerical p8| , p9| , p3| , p5[ and analytical pC| . 

In Sect. Ill we have developed an effective continuum 
field theory for the Haldane phase. The results obtained 
there can be commented as follows: 

i) At T = we obtain a spin gap and a vanishing 
magnetization up to the lower critical field H Cl . At the 
mean-field level the elementary excitation spectrum con- 
sists in a triplet of magnons, and one magnon branch 
goes soft at H = H Cl , thus signalling an instability and 
the transition to a gapless regime. Parenthetically, in a 
recent paper Loss and Normand |^j have developed an 
approach quite similar to ours. They claim however that 
the magnetization does not vanish in the spin-gap phase 
at the mean-field level. On the light of the results of 
Sect. Ill such a claim seems to be incorrect, and a com- 
pletely consistent picture for the Haldane phase appears 
to emerge already at the simplest mean-field level. In its 
present form and with an appropriate redefinition of the 
coupling constant and of the spin- wave velocity |i"o| , the 
theory could apply also to the description of spin lad- 
ders in an external field [ |37| in the regime in which the 
correlation lenght is much larger than the width of the 
ladder 

ii) It is interesting to remark that no singularities can 
develop at finite temperatures. Indeed, as already re- 



marked immediately after Eq.([75|), setting £~ 2 = the 
last term in (^) would develop a severe infrared diver- 
gence for any (3 < oo (strictly) which is altogether absent 
at T = 0. Therefore at T > £ remains finite and the 
system is in a disordered (gapped) phase, for all values 
of the applied field, as it should. To be more precise, 
finiteness of £ implies A > H (strictly) for all values of 
H, with the following consequences: a) no magnon mode 
will go soft at any T > as, in that case (see Eqs.([77|) 
and (|S0|)), the spectrum of elementary excitations will be 
given by: e = eo,e± with: e± = y A 2 + c 2 q 2 ± H and 
hence: e_(g = 0) = A — H > (although exponentially 
small), and: b)the magnetization (see Eq. (p7[) ) remains 
smooth all the way up to saturation, although exponen- 
tially small at low temperatures. 

What we have presented in Sect. Ill is a field-theoretic 
approach possibly at its simplest level. At least for fields 
H < H Cl , however, the role of quantum fluctuations can- 
not be underestimated and should be taken into account. 
On this direction, as well as in that of obtaining a consis- 
tent field-theoretic description of the gapless phase, work 
is presently in progress and we hope to report on it in 
the near future. 
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